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We study maximally multipartite entangled states in the context of Gaussian continuous variable quantum 
systems. By considering multimode Gaussian states with constrained energy, we show that perfect maximally 
multipartite entangled states, which exhibit the maximum amount of bipartite entanglement for all bipartitions, 
only exist for systems containing n = 2 or 3 modes. We further numerically investigate the structure of these 
states and their frustration for n < 7. 
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I. INTRODUCTION 

Entanglement is nowadays recognized as a fundamental re- 
source for quantum information processing (see e.g. |T]). It 
explicitly appeared long before the dawn of quantum infor- 
mation science and without any reference to discrete variables 
(qubits) |2]. In fact, it first came to light in the context of 
continuous variables [3]. Thus, its characterization must nec- 
essarily include the latter as well. Along this line, important 
milestones have appeared in terms of continuous variables and 
more specifically Gaussian states (see e.g. 0] and reference 
therein). 

Although bipartite entanglement can be conveniently char- 
acterized (e.g. in terms of purity or von Neumann entropy) fsl, 
the characterization of multipartite entanglement remains a 
challenging problem, together with the definition of a class of 
quantum states that exhibit high values of multipartite entan- 
glement. Recently, the notion of maximally multipartite en- 
tangled state (MMES) was introduced in the qubit framework 
10]. These states have a large (in fact, maximum) value of av- 
erage bipartite entanglement over all balanced bipartitions of 
a system of qubits |6, 7]. They are solution of an optimiza- 
tion problem and minimize a suitably defined cost function, 
that can be viewed as a potential of multipartite entanglement. 
A MMES is called "perfect" if its average entanglement satu- 
rates the maximum bipartite entanglement for all bipartitions. 
Perfect MMESs exist for n = 2, 3, 5 and 6 qubits, they do not 
exist for ri = 4, n>7||^|8t], while the case n = 7 is still an 
open problem. In terms of potential applications in quantum 
information science, MMESs are the ideal resource for ini- 
tializing a quantum internet and could be useful in several 
multiparty quantum information protocols (like e.g. controlled 
teleportation fiol] or quantum secret sharing 111]). 

The concept of MMES was extended to the framework of 
continuous variable (and Gaussian) systems in (l^. There a 
Gaussian MMES is a state with a maximal rank for any bipar- 
tition of the n party system in the limit of infinite squeezing 
fl^. Notice that such a state allows perfect quantum tele- 
portation among its n parties. Here, with the aim of char- 
acterizing Gaussian MMESs, we adopt a different viewpoint, 
by introducing a constraint on the maximal mean energy al- 



lowed per user (which, eventually, will be let to go to in- 
finity). Hence, we look for states that present the maximal 
amount of bipartite entanglement compatible with the given 
constraint. We will show that, following this definition, per- 
fect MMESs exist in the continuous variable Gaussian set- 
ting only for n — 2, 3. For n > 4 a simple argument 
shows that perfect MMESs do not exist, hence manifesting the 
phenomenon of entanglement frustration [13] (see also [14]). 
Then, for 4 < n < 7 we study the distribution of entangle- 
ment among the bipartitions. Finally, we find examples of 
MMESs and provide numerical evidence that bipartite entan- 
glement can be optimally distributed for n = 5,6. 

II. BASIC DEFINITIONS 

A system composed of n identical (but distinguishable) 
subsystems is described by a Hilbert space 7i = Tis, with 
T^s '■= ®ies '^^'^ ^ ~ {1; 2, ... , n}, which is the tensor 
product of the Hilbert spaces of its elements t)i ~ f). Exam- 
ples range from qubits, where (j = C^, to continuous variables 
systems, where f) — i^(K). We will denote a bipartition of 
system S by the pair (A, A), where AcS, A = S\ A and 
1 < JT-A < 'with UA — \A\, the cardinality of party A (of 
course, ua + = n). At the level of Hilbert spaces we get 

n^HA^HA- (1) 

A crucial question in quantum information is about the 
amount of entanglement between party A and party A. When 
the total system is in a pure state E H, which is the only 
case we will consider henceforth, the answer is simple, and 
can be given, for example, in terms of the purity 

TTA = tr(pi) (2) 

of the reduced density matrix of party A, 

PA:^tTnA\^)m/M'. (3) 

Indeed, this quantity can be taken as a measure of the entan- 
glement of the bipartition {A, A). Its range is 

<A <^A<1, (4) 



2 



where 



<nt. = (diniWA)-^ = (dim())- 



(5) 



with the stipulation that l/oo = 0. 

The upper bound 1 is attained by unentangled, factorized 
states IV') = \4>)a ® Ix)a- On the other hand, when dimi) < 
oo, the lower bound, which depends only on the number of 
elements ua composing party A, is attained by maximally 
bipartite entangled states, whose reduced density matrix is a 
completely mixed state 



PA 



(6) 



Note, however, that for continuous variables, dim \) = oo, 
the lower bound t^^^ = is not attained by any state. There- 
fore, strictly speaking, in this situation there do not exist max- 
imally bipartite entangled states, but only states that approx- 
imate them. This inconvenience can be overcome by intro- 
ducing physical constraints related to the limited amount of 
resources that one has in real life. This reduces the set of 
possible states and induces one to reformulate the question 
in the form: what are the physical minimizers of (|2|i, namely 
the states that minimize (O and belong to the set C of physi- 
cally constrained states? In sensible situations, e.g. when one 
considers states with bounded energy and bounded number of 
particles, the purity lower bound 



inf{^^, IV-) e C} > ^; 



riA 
min 



(7) 



is no longer zero and is attained by a class of minimizers, 
namely the maximally bipartite entangled states. If this is the 
case, we can also consider multipartite entanglement and ask 
whether there exist states in C that are maximally entangled 
for every bipartition {A, A), and therefore satisfy the extremal 
property 



nA,C 



(8) 



for every subsystem ^ C S'withn^ = < n/2. In analogy 
with the discrete variable situation, where dim fi < cx) and 
C = we will call a state that satisfies ^ a perfect MMES 
(subject to the constraint C). 

Since the requirement ([8]) is very strong, the answer to the 
quest can be negative for n > 2 (when n = 2 it is trivially sat- 
isfied) and the set of perfect MMES can be empty. We remind 
that for a system of n qubits, i.e. when dim t) = 2, perfect 
MMESs exist for n = 2, 3, 5, 6, do not exist for n = 4, n > 7 
|6, 8], while the case n = 7 is still open. This is a symptom of 
frustration 1.13.1 . We emphasize that this frustration is a con- 
sequence of the conflicting requirements that entanglement be 
maximal for all possible bipartitions of the system. 

In the best of all possible worlds one can still seek for 
the (nonempty) class of states that better approximate perfect 
MMESs, that is states with minimal average purity. We there- 
fore consider the potential of multipartite entanglement 



ttme = ^{t^a) ■= 



n 
[n/2] 



\A\ = [n/2] 



TTA 



(9) 



where the sum runs over all balanced bipartition, with \A\ = 
iiA — [n/2], where [ • ] denotes the integer part. (It is immedi- 
ate to see that a necessary and sufficient condition for a state 
to be a perfect MMES is to be maximally entangled with re- 
spect to balanced bipartitions, i.e. those with ua — [?^/2].) By 
definition a MMES is a state that belongs to C and minimizes 
the potential of multipartite entanglement. Obviously, when 

""me. min ■— mine TTME = ""rnhf' ''^ there is no frustration and 
the MMESs are perfect. In order to quantify the amount of the 
frustration (for states belonging to the set C) we will take the 
quantity 



Xn 



'ME.i 



Jtt 



n/2],C 
mill 



(10) 



Eventually we will consider the limit C ^ Ti,. 



III. BASIC TOOLS FOR GAUSSIAN STATES 

Let us consider a collection of n identical bosonic oscilla- 
tors with (dimensionless) canonical variables {qk,Pk}k=i....n- 
We assume that the oscillators have unit frequency and set 
^ = 1. A quantum state of the n oscillators can be described 
by a density operator on the rt-mode Hilbert space, or 
equivalently by the Wigner function on the rt-mode phase 
space 

w^(n) (q, P) - y d"y(q - y|P(„) lq + y)e'"^-p, (11) 



where q := (91,... g„), p := {pi, 
{yi, . . .i/n) G K", and we have denoted by 



|q±y> 



ik=i\qk±yk) 



■Pn), y := 



(12) 



the generalized eigenstates of the 'position' operators qk- By 
definition, Gaussian states are those described by a Gaussian 
Wigner function. Introducing the phase-space coordinate vec- 
tor X = {Xi,...X2n) ■= (gi, Pi, ... a Gaussian 
state has a Wigner function of the following form: 



W^(n)(X) 



1 



(27r)"^det(V) 
X exp 



-i(X-Xo)V-i(X-Xo)T 



,(13) 



where Xq = (X), with (/(X)) / /(X)W^(„) (X)d2"X, is 
the vector of first moments, and V is the 2n x 2n covariance 
matrix (CM), whose elements are 



Yl,n = {{Xi - {Xi)){X„, - {X„,))). 



(14) 



We will also consider an equivalent representation de- 
fined by a different ordering of the canonical variables X = 
{qi,q2, ■ ■ ■qn,Pi,P2 ■ ■ -Pn)- In this representation the CM is 
denoted V and has elements 



^ {{Xl - {Xi))iXrn - {X^))). 



(15) 
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In order to study the properties of entanglement for Gaus- 
sian states we will consider the purity 



'^(P(n)) tr(p^„)). 



(16) 



From Eq. ( fTTT l. it is straightforward to compute this quantity 
in terms of the Wigner function: 

^(P(„)) - (27r)" / [M^(„)(X)]'d2nx. (17) 



In particular, the purity of Gaussian states is a function of the 
determinant of the CM. From Eqs. ( fTTI l it follows that 



1 



2'Vdet(V) 



(18) 



with the bound Tr{p(^n) ) < 1- Notice that from ( fTSI l a Gaussian 
state with positive CM is pure if and only if 



The inequality ( l23T l can be written in terms of the vector of 
first moments and the trace of the CM as follows 



tr(V) + |(X)|- 
2n 



<N+l. 



(24) 



Now we notice from Eq. (fTSl l that the Gaussian state mini- 
mizing the purity is the one whose CM has maximal determi- 
nant under the constraint. The determinant and the trace of a 
CM are functions of its eigenvalues {wj}j=i,...2n- We hence 
consider the problem of finding V such that det (V) = Y\j ""j 
is maximal under the constraint tr(V) = ''^j ^ 2n{N + 
1/2) — |(X)p. The unique solution is the scalar matrix [N + 
1/2 — |(X}|^/(2n)]l2„. For a given value of (X) the maximal 
value of the determinant is hence [A^+ 1/2 - | (X) |2/(2n)]2". 
It follows that the unique Gaussian state minimizing the purity 
under the energy constraint ( |23] ) — and hence the constraint 
(I20I 1 — is the one with (X) = and CM as in Eq. (|2T]i, i-e. 
the thermal state. QED 



detV= ( - 



(19) 



IV. GAUSSIAN MMES 



As anticipated in Sec. |ll] in order to obtain sensible results, 
we will impose impose a suitable energy constraint. Here we 
do not allow more than N mean excitations for each bosonic 
mode, i.e. 

ilk+Pk) < 1 foj. k^l,...n. (20) 
2 ~ 2 

This constraint introduces a cutoff in the Hilbert space of each 
quantum oscillator 

A particular example of Gaussian state is the thermal state 
p^^^ with N thermal excitations per mode described by a 
Gaussian Wigner function with vanishing first moments and 
CM 



V"^ = (iV + 1/2)I2„ 



(21) 



Obviously, p*^^ satisfies the constraint ( l20l i. We now show the 
following 

Proposition 1 Among all Gaussian states, the thermal state 
is the unique state that minimizes purity under the constraint 
The corresponding minimal purity is 



1 



{n)l 2"(iV+l/2)"' 



(22) 



Proof: We prove that the thermal state is the unique mini- 
mizer of the purity among the Gaussian states satisfying the 
inequality 



1 vM±zl) 

2 



(23) 



which constrains the average mean energy per mode. This is 
indeed sufficient to prove the proposition since all the states 
satisfying (l20l i also satisfy ( l23T l. 



In the following, we will focus our attention on the case of 
pure Gaussian states characterized by Eq. (fT9] ) and subjected 
to the energy constraint (|20] |. As in Sec. Ull we consider a 
collection of n-modes and a bipartition into two disjoint sub- 
sets A and A, containing ha and riA — n ~ ha modes, re- 
spectively, with HA < riA- In order to quantify the bipartite 
entanglement between the two subsets of oscillators, we com- 
pute the purity tta of the reduced state of subsystem A. The 
modes of subset A and A have phase-space coordinates 'Ka 
and X^, respectively. The Wigner function describing the re- 
duced state of party A is obtained by integrating the Wigner 
function of the whole system over the variables belonging to 
A, i.e. 



Ty(„)(x) d'^-l^A■ 



(25) 



It follows that the reduced state of a Gaussian state with first 
moment (X) and CM V is Gaussian with first moment (X^) 
and CM Wa- The CM of the reduced state is the square sub- 
matrix of V identified by the indices belonging to subsystem 
A. 

States that are maximally entangled with respect to the 
given bipartition are those admitting a reduced state for sub- 
system A with minimum value for the purity. From Proposi- 
tion [T] the reduced system has to be in a thermal state of ua 
oscillators. Taking into account the energy constraint ( |20| ) we 
get 



\ \tp) en, Gaussian, 



{qI+pD 



1 



2 - 2 ' ' 



and 



TT • ^ 7r4 ^ i 
mm — ^ — 



(26) 
(27) 



•4.U "^A - C nA,N • -1 

with7r^i„ =7r,^i^ given by 
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We will generalize this property to multipartite entangle- 
ment by requiring minimal possible purity for each subsys- 
tem A of the modes — assuming that the state of the total 
system is pure — thus defining a Gaussian maximally multi- 
partite entangled state. In particular, we will define a perfect 
Gaussian MMES as a pure Gaussian state of n oscillators that 
is maximally entangled with respect to all balanced biparti- 
tions and satisfies the energy constraint ( |26] ). It follows from 
this definition that the reduced state is thermal for all possible 
bipartitions (^,^4). 

In order to formalize the above definition of MMES for 
Gaussian states with constrained energy we start from (|9]l and 
define a (normalized) potential of multipartite entanglement: 



X ttme/tt. 



mill 



(Ar + l/2)"-^E dct(VA) 



\-i/2 



(28) 



where ua — [n/l], N a is the square sub-matrix defined by the 
corresponding indexes. The minimum of this quantity 
(where C ~ is the constraint) is a measure of frustration, 
according to Eq. dTol l. A Gaussian MMES will be a minimizer 
of the potential ( |28] ). The potential in ( |28] ) is the normalized 
purity of the reduced state, averaged over all balanced biparti- 
tions. Notice that X > 1, and perfect Gaussian MMESs satisfy 

We recall from the above discussion that the requirement 
of the minimization of purity for a given bipartition could be 
in contrast with that for another bipartition. Therefore perfect 
MMESs do not necessarily exist. Actually, we get the follow- 
ing 

Theorem 1 Perfect Gaussian MMESs only exist for n < A. 

Proof. First we present examples of Gaussian MMESs for 
n = 2,3, then we show that Gaussian MMESs do not exist 
for n > 4. 

A two-mode Gaussian state is described by the vector of 
first moments (X) = (((71,^1,(72,^2)) and by the 4 x 4 CM 



Vi,i 



Vl,2 

¥2,2 



(29) 



Imposing that the one-mode reduced states are thermal im- 
plies (X) = (0, 0, 0, 0) and Vij = ¥3,2 = [N + l/2)l2. It 
remains to specify the submatrix Vi_2 in order to obtain a well 
defined CM V satisfying the purity condition ( fT9] l. A solution 
is given by the CM 



/ cosh r sinh r 

cosh r — sinh r 

sinhr coshr 

\ — sinh r cosh r 



(30) 



describing a two-mode squeezed state, the so-called twin- 
beam state 101, for 



coshr = 2A^+ 1 



(31) 



(squeezing parameter r/2). 

Let us now consider the case n = 3. A three-mode Gaus- 
sian MMES has all the three single-mode reduced systems in 



a thermal state. Hence the vector of first moments vanishes, 
and the CM has the form 

[N 



+ l/2)l2 

VI2 

vT,3 



V1.2 V1.3 

+ 1/2)I2 V2.3 

(iV + l/2)l2 



(32) 



It remains to determine the sub-matrices Vi 2, Vi^a, ¥2,3 in 
order to obtain a well defined CM obeying the purity con- 
strain ( fT9] l. A solution is given by an instance of the tripar- 
tite Gaussian GHZ states 1II6II characterized by the condition 

Vi,2 = Vi^a = ¥2^3 = diag(z;+, w_) with 



v± 



N{N +1) 
AN + 2 



1± 



(4iV + 2)2 
2N{N+l) 



(33) 



The cases n — 2, 3 are the only ones in which perfect Gaus- 
sian MMESs exist. The non existence of perfect Gaussian 
MMESs for n > 4 is easily seen by inspecting the n-mode 
CM. Indeed, the generic submatrix of n^i < n/2 modes is of 
the form 



V' 



(34) 



J. For 



where is a 2n^ x 2nA matrix and A = {ii, . . . 
n > 4, the definition of perfect Gaussian MMES implies that 

Vij = for i,jeA and i ^ j, (35) 
V,,, = {N+ l/2)l2 for i e A. (36) 

This condition must hold for all bipartitions {A, A) and, there- 
fore, all off-diagonal sub-matrices are zero. As a conse- 
quence, the CM of the Gaussian state is diagonal of the form 
(N + 1/2)I2„. Such a CM describes a thermal state with N 
thermal excitation per mode, in contradiction with the require- 
ment that the global state of the n oscillators is pure. QED 



V. NUMERICAL SEARCH OF MMES 

We have seen that perfect MMESs only exist for n < 4. 
For n > 4 we now numerically search for 71-mode pure states 
minimizing the cost function ( |28] ), under the energy constraint 
( I20I ). Minimizing the cost function ( |28] l corresponds to min- 
imizing the average purity of the reduced states. The value 
X = 1 corresponds to a perfect MMES. For n > 4 this is 
possible only for iV = 0, where the problem becomes trivial 
since the only state compatible with the energy constraint is 
the vacuum, which is a separable state. 

For numerical investigations we use a convenient 
parametrization of n-mode pure Gaussian states. First 
of all in the following we will assume, without loss of gener- 
ality, (X) = (X) — 0. It remains to parameterize the set of 
n-mode covariance matrices. Working in the representation 
( fT5l ), it is possible to show that the CM of a n-mode pure 
Gaussian state can be written as ITtIi 



1, 



IT 



2mT 



(37) 
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FIG. 1: Minimal value of the (dimensionless) cost function x ver- 
sus the (dimensionless) energy constraint parameter A'^, for several 
values of n. 



where T is a diagonal matrix of the form 

; 



T = 







(38) 



with K a non singular diagonal matrix, while K is both sym- 
plectic and orthogonal. Therefore it has the form 



X Y 
-Y X / ' 



(39) 



where the matrix U = X + iY is unitary. 

Figure [T] shows minimal value Xmin °f '^^e potential 
under the constraint ( |20l ). for 4 < n < 7, as a function of 
the mean number of excitations per mode N. This minimal 
value yields a measure of \he frustration present in the system, 
which does not allow the existence of perfect MMESs. The 
larger the minimal value of x, the larger the frustration. The 
numerical analysis indicates that the minimum of the poten- 
tial of multipartite entanglement is a nondecreasing concave 
function of N; moreover, a plateau is reached for sufficiently 
high values of A^. This saturation value increases with n, but 
oscillates between even and odd n. 

Since it is not possible to find perfect MMESs, it is im- 
portant to quantify the distribution of entanglement. A good 
distribution of entanglement should be rather insensitive to a 
change of bipartition. A fairly distributed multipartite entan- 
glement should therefore be characterized by a distribution 
(over balanced bipartition) with a small standard deviation 
111811 . We therefore consider the standard deviation of the pu- 
rity of the reduced states over balanced bipartitions 



Ax ■.= ^J{N + 1/2)^"-^ E [det(VA)-i] - 



(40) 



with ua = [n/2]. We will call a MMES with Ax = a uni- 
formly optimal MMES, because it has an optimal distribution 



FIG. 2: Standard deviation Ax (dimensionless) of the normalized 
purity versus the (dimensionless) energy constraint parameter A'^, for 
n = 4 and 7. For n — 5 and 6 the MMESs have optimally distributed 
entanglement, hence vanishing standard deviation. 



TABLE I: Comparison between qubit and Gaussian maximally mul- 
tipartite entangled states. 



n 


qubit perfect MMES 


Gaussian perfect MMES 


2,3 


yes 


yes 


4 


no 


no 


5,6 


yes 


no, but uniformly optimal* 


7 


no 


no 


> 8 


no 


no 



'numerical evidence 



of entanglement: entanglement (and frustration) is fairly dis- 
tributed over all bipartitions that attain the minimal value of 
purity allowed by frustration. Of course, a perfect MMES is 
uniformly optimal. 

Figure |2]displays the pecuUar behavior of the standard de- 
viation of the purity: such standard deviation has a different 
behavior as a function of N for different values of n. For 
n = 5 and 6, MMESs are not perfect. Nonetheless, inter- 
estingly enough, they have an optimal distribution of entan- 
glement: Ax = (for n = 5,6). By contrast, a non opti- 
mal entanglement distribution has been found for MMES with 
n = 4 and 7. This reminds, mutatis mutandis, of what hap- 
pens in qubit systems, where frustration appears for n = 4 
and n > 7 l^llllMl- The numerical analysis shows a Ax of 
the non uniformly optimal states which is a concave nonde- 
creasing function of the energy parameter N. We notice also 
in this case the presence of a saturation effect for large values 
of N. These findings are summarized in Table U 
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VI. CONCLUSIONS 

In conclusion, we have characterized Gaussian states that 
display a maximal amount of multipartite entanglement com- 
patible with a given constraint on the mean energy. We have 
shown that perfect Gaussian MMESs (that saturate the max- 
imum mean energy) only exist for n = 2 and 3, while the 
phenomenon of entanglement frustration appears already for 
n> A. Curiously, we found clear numerical evidence that al- 
though perfect Gaussian MMESs do not exist for n = 5 and 
6, for these particular values of n, bipartite entanglement can 
be optimally distributed, in the sense that the standard devia- 
tion of purity over balanced bipartitions (|40] | can be made to 
vanish. We numerically found that, by contrast, such standard 
deviation cannot be made to vanish (and bipartite entangle- 
ment is therefore not optimally distributed) for n = 4 and 7. 
This peculiar situation is reminiscent of that encountered with 
qubit MMES, where perfect MMESs exist for n = 5 and 6 
(qubits), but do not exist for n — A and (probably lj6|]) n = 7. 
This suggests once more that n — 2,3,5 and 6 are "special" 
integers. We endeavored to summarize these conclusions in 
TableU] Experience with integers does not induce us to expect 
that these amusing peculiarities only occur for n < 6 (for in- 
stance, we numerically found that uniformly optimal MMESs 
— with vanishing Ax — exist for n = 9). Additional research 



is needed in order to investigate the large n behavior and clar- 
ify the underlying structure of entanglement frustration. We 
emphasize again that this frustration is a consequence of the 
conflicting requirements that entanglement be maximal for all 
possible bipartitions of the system. The same phenomenon is 
also worth studying under different constraints, like, for in- 
stance, the (weaker) energy constraint ( |23] l. 

From a more applicative perspective, we emphasize that 
due to recent progress in the optical generation of Gaussian 
entangled states (up to 9 modes |19]) the above features are 
also liable to experimental check. These results, combined 
with those obtained in lfl2ll . and the ensuing proposed char- 
acterization of entanglement would help in optimizing mul- 
tiparty quantum information protocols with continuous vari- 
ables. 
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